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Abstract. - We have applied the recently developed dual fermion technique to the spectral prop- 
erties of single-band Anderson impurity problem (SI AM). In our approach a series expansion is 
constructed in vertices of the corresponding atomic Hamiltonian problem. This expansion con- 
tains a small parameter in two limiting cases: in the weak coupling case (U/t — > 0), due to the 
smallness of the irreducible vertices, and near the atomic limit (U/t — > oo), when bare propagators 
are small. Reasonable results are obtained also for the most interesting case of strong correlations 
(U ~ t). The atomic problem of the Anderson impurity model has a degenerate ground state, so 
the application of the perturbation theory is not straightforward. We construct a special approach 
dealing with symmetry-broken ground state of the renormalized atomic problem. Formulae for 
the first-order dual diagram correction are obtained analytically in the real-time domain. Most 
of the Kondo-physics is reproduced: logarithmic contributions to the self energy arise, Kondo-like 
peak at the Fermi level appears, and the Friedel sum rule is fulfilled. Our approach describes also 
renormalization of atomic resonances due to hybridization with a conduction band. A generaliza- 
tion of the proposed scheme to a multi-orbital case can be important for the realistic description 
of correlated solids. 



Introduction. — The problem of realistic description 
of spectral properties of correlated impurity in the metal- 
■ lie surrounding is far from solution despite of considerable 
progress during the past years. There are quite success- 
ful tools for large Coulomb U , particularly diagrammatic 
pseudoparticle approach [1]. Transport properties can be 
described in this limit also [2,3]. The case of moderate 
U is less studied, although there is an interesting physics 
behind, e.g. renormalization of the multiplet structure 
due to hybridization with conduction band. An exact 
solution by using the Bethe-ansatz [4, 5] is obtained for 
thermodynamical properties only. There are numerically 
exact continuous-time quantum Monte Carlo calculations 
in imaginary time domain [6,7]. Practical accuracy of 
these methods is limited by an ill-posed problem of ana- 
lytical continuation to the real-frequency axis. More ac- 
curate data can be obtained within the Numerical Renor- 
malization Group framework [8]. Finally, it is important 
to mention a simple RPA-like approach which starts from 
the local-moment broken symmetry solution [9] and de- 
scribe Kondo properties surprisingly well. Abrikosov-Suhl 



resonanse has been reproduced properly, it fulfils Friedel 
sum rule and its width corresponds to the Kondo energy 
scale. 

In the present paper we describe an approach aiming 
analytical description of spectral properties of correlated 
impurity with moderate U. In this first paper we apply 
our method to a single-band Anderson impurity model 
(SIAM). Our approach, which starts from the exact map- 
ping of SIAM to auxiliary (dual) variables [10,11], resem- 
bles key properties of the broken-symmetry self-consistent 
RPA approach [9], but allows fully analytical consider- 
ation. Another important property of our method is a 
renormalization procedure, which allows to describe a shift 
of the atomic resonances due to hybridization with a con- 
duction band. 



Anderson impurity problem. We consider a 
single-band Anderson impurity problem at zero temper- 
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ature. The action is 



I- /-t-OC 

S = S at - dt dt'c at A(t - t')c at > 

J J — oo 

S at = J dt ^ic at ^c at - Un^ t n lt + ^h a ^j 

Summations over spin index a =t,| are supposed; fi at = 
^(cat-aCat + Cat+oCat) corresponds to n — \ term in the 
Hamiltonian [12], so that there is a particle- hole symmetry 
of Sat if chemical potential n equals zero. 

To build up the theory, we suppose an adiabatic switch 
of the hybridization at infinity. Formally this means that 
A is multiplied by a slow profile function A prof (t) such as 
AP rof (±oo) = and AP rof (t) = 1 otherwise. So, the evo- 
lution starts and finishes with a bare atom having a single 
electron. The complete evolution operator S(— oo, oo) pre- 
serves the spin orientation. Therefore, it can be divided 
into the two parts, responsible for the evolution that starts 
and finishes with a certain spin orientation: 



S = S 



rr 



s 



u. 



(1) 



where iS-f-f andSj^ describe the evolution | T)-oo — » | T)+oc 
and | l)-oo — » | 4)+oo- Formally, S a a can be defined as 



rW) 

J\a) 



V[cc], 



(2) 



where V[cc] implies integration over the trajectories 
starting and finishing with a definite spin orientation. 
Such a trick is necessary because of the degeneracy of the 
ground state of the atomic Hamiltonian. Remind that 
for usual stationary diagrammatic technique, the ground 
state is unique, and therefore always evolves to itself at 
infinite time. 

In practice, it is enough to consider only one part of 
the evolution operator, for example S^. Clearly, the ac- 
count of another part is equivalent to spin-averaging of the 
results obtained. 

Preliminary analysis. — A transformation to the 
dual variables [11] requires a splitting of the action into 
two parts. The first part might be nonlinear but should 
be exactly solvable, and the second part should be Gaus- 
sian. The simplest way is to choose S a t and hybridization 
as those parts. Such a theory is described in the next two 
sections. It will be shown that it describes the low-energy 
physics quite good. However, it will be shown also that 
a correct description at all frequencies requires a more 
sophisticated approach, and a renormalization procedure 
will be introduced. 

Hubbard-Stratonovich decoupling of the hybridization 
and subsequent integration over c, c in the evolution oper- 
ator results in the dual action 



S[f,f] 



+ 00 



9**, e )f*,cU,* + V[f, /], 
(3) 





Fig. 1: The world-lines contributing to g aa (a,b) and 7^ Ti (c). 

where Taylor coefficients of the nonlinear part V[f,f] are 
vertexes of the atomic problem. 

Momenta of the atomic problem can be calculated in 
the time domain, using world-line representation. Parts 
(a) and (b) of Figure 1 show the nonvanishing world lines 
describing the calculation of the Green's function. One 
can see, that since the evolution starts and finishes with 
a single electron with the spin-up orientation, there is 
a definite ordering of the creation-annihilation operators: 
annihilation should precede creation for spin-up and vice 
versa for spin-down operators. At half-filling one obtains 
g n (t) = i6{-t)e i( - u / 2 » and g u (t) = -idifye^l^ 1 (0 is 
a Heaviside step function). Fourier transform gives 



5TT ( e ) = 



(4) 



e + U/2-iO' ' e-U/2 + iO' 

These formulas are sufficient to construct the mean-field 
approximation, as higher momenta docs not enter the the- 
ory in that case. Very similarly to formulas (14-15) of [11], 
one obtains 



g^\e)=G IJU {e)-g aa {e) 



(5) 



(through the paper, the mean-field results are denoted 
with the calligraphic letters). 

Let us remind that an average over spin indices should 
be taken to obtain the final result, so we get 

1/2 1/2 



+ U/2 - A(e) - iO e - U/2 - A(e) + iO 



(6) 



(In our notation, quantities before and after spin averaging 
are marked with the same letters. To escape confusion, we 
always supply non- averaged quantities with spin indices.) 

We will be mostly interested in Kondo-likc problem, so 
that atomic resonances ±17/2 lie outside the conduction 
band. One can see the mean-field DOS of such a problem 
is built of slightly reshaped band DOS and two i5-peaks at 
±{7/2. The main fault of this approximation is that no 
Kondo peak appears near Fermi level, while the higher- 
energy part of DOS is qualitatively correct. 
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Low-energy properties: general consideration. 

A calulation of the corrections to the mean-field theory 
requires knowledge of higher momenta of the atomic prob- 
lem. The two-particle Green's function 01234 = (C1C2C3C4) 
can be calculated similarly to 512 (indices here are com- 
binations of energy and spin, for instance 1 stands for 
ei, <7i). The expression for contains 24 different terms, 
corresponding to various mutual order of the four time 
and spin arguments. However it turns out that only 
four of whose terms contribute to the non-Gaussian part 

ri234 = 0I234 - 013524 + 014023- The world lines corre- 
sponding to these terms are depicted in part (c) of Fig- 
ure 1. After the Fourier transform (omiting the energy- 
conserving delta function), we obtain a simple formula 

for the fourth-order vertex 7^34 = — *0i~i 1 022 1 033 1 044 1 -^i234- 
Four-point vertex for all-the-same spin indices equals zero, 
and for the different indices 



7{tu( ei ' e2 ' e3 ' e4 ) = ~ U 



U 2 



£3 - e2 — Oi 



(7) 



The first term of this expression is local in time, whereas 
the second one is proportional to 6^3 — ti)5(t4 — ti)S(t 3 — 
£2) in time-domain. The Heaviside function appears here 
because of the degeneracy of the ground state: a time 
interval between the pairs t\ti and £2^3 can be arbitrary 
large (see Figure 1). 

In this paper, we restrict ourselves with the simplest 
approximation beyond mean-field, that is the first-order 
diagram correction to the dual self-energy: 



^dual 



2tt 



(4) 1 



• £ > e Jt/J4 ( e > de 



(8) 



and similarly for Sf" a '. It would be more accurate to use 
the integrand with the renormalized dual Green's function 
Gf^ al , but we will stay with the expression with Gfl al for 
simplicity. 

To obtain an expression for the Green's function, it's 
practical to use the identity (18) of [11] . We obtain 

G~(e) = — —± — - 

g^(e) A(e) - (g aa {e) + (S^'(e))- 1 ) 

(9) 

An important property of these formulas is 
that S^ ai contains the Kondo logarithm. In- 
deed, substituting (7) into (8) and taking into ac- 
count that lmQ dual changes its sign at the Fermi 
level, one finds that Re£^" a/ contains the logarith- 
mic singuliarity: (27T)- 1 J e _^_ i0 Imgffi ai (e')<fe' « 
7r- 1 J7 2 log(-O/e)Im^J ta/ (-0), where the cut-off ft is 
about the half-bandwidth. Futher substitution of Eq. (5) 
gives an estimation ReS^" ai w 2n~ 1 log(— Q/e) for the 
case of large U . Clearly, such a logarithmic behaviour is 
also reflected in DOS. 

The dual perturbation theory is valid while the correc- 
tions to the mean-field result are small enough. Accord- 
ing to formula (9), it means that the domain of validity 



is determined by the inequality 0-fT^Tt ^ 1- Note that 
at the Kondo energy [13] ck cx p(~ ^) the left-hand 
side of this inequality equals one-half, so our theory is 
formally valid only above the Kondo scale. But it is im- 
portant to note that some of low-energy physics is also 
reproduced. Namely, Fricdel sum rule [13] is satisfied: as 
one can see from (9), the divergence of S^ a/ corresponds 
to G aa = -A" 1 . 

A calculation for a concrete system should be presented 
to ensure that the theory behaves well. In the next sec- 
tion, we present the analytical formulas and corresponding 
graphs for the case of half-filled system in a semicircular 
bath. 

Low-energy properties: analytical results for a 
semicircular bath. — An explicit calculation of the 
first-order diagram in dual variables is a simple although 
not trivial task. The bare dual Green's function Q^ al (e) 
has an atomic pole at — U/2 (and at U/2 for Gff (c)) as 
well as a branch cut originating from the fact that A(e) 
represent a system with a continuous spectrum. 

In our calculation we will use a semicircular A(e): 



A(e) = 



1 (2t 



2 \D 



(e — sgn(e)\/ e 2 — D 2 ) 



(10) 



Here D is the half-bandwidth, which for simplicity is put 
to one in all calculations below. Hybridization constant t 
describes the coupling between the impurity atom and its 
nearest neighbor. 

The diagram we are going to calculate is an integral over 
the real axis: 



u 2 



e-e'-iO 



Gfl al {e') (11) 



de' 



U 2 



Efr'to = / ^ ( u + e ,_ e _, ) sfrV) (12) 



Let's examine the pole structure of the part proportional 

due 

rr 



to U (constant part) in the expression for E^ a/ (e) . Denote 



this part S*" 2 ' 
U 



Sdual 
u — ~ l 



2TT 



Iv 



A(e') de 1 



U/2 + i0)(e' - U/2 + i0- A(e')) 

(13) 

There are a pole at e' = U/2 — iO inherited from the 
atomic problem and at most two poles which are solutions 
of the equation e — U/2 = A(e). For our special choice of 
A(e) this equation is reducible to the following system: 



ReA(e) 



U/2 



ImA(e) =0 

e 2 (l - At 2 ) + eU(2t 2 - 1) + ((2t 2 ) 2 + (U/2) 2 ) = 
|e|>l 

(14) 

In the case of t = 1/2 this quadratic equation turns 
to a linear one and its root eo = \ (U + y) is always 
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Fig. 2: Graphical solution of the equation Re A(e) = e — U/2. 

greater than or equal to 1. When t is not equal to 1/2, two 
real roots e + and e_ are possible. This is a manifestation 
of an additional energy level splitting introduced by the 
"defect" of the lattice. However in the limit U 3> 2t the 
only relevant root is still e+. To understand this statement 
it's useful to represent ReA(e) = e — U/2 equation in a 
graphical form (see Figure 2). 

As one can see, the magnitude of t determines the height 
of "shoulders" in this figure, while U/2 is a vertical dis- 
placement of the diagonal line e— U/2. At U large enough 
in comparsion with 2t 2 (dimensionless) there is only one 
intersection of the curves on the positive half-axis, i.e. e + . 
A complete analisys of the equation yields the following 
results: 



t,u 


roots 


t = 1/2 

t ^ 1/2, 2t 2 < l,U/2 < 1 - 2t 2 
t ^ 1/2, U/2 > |1 - 2t 2 \ 
2t 2 > 1, U/2 <2t 2 - I 


no real roots 

6 = 6+ 

e = £+,£_ 



{U/2){U 2 -2)± 4tV4t 2 + {U/2f 
2(4t 2 - 1) 



(15) 



Both poles U/2 — iO and e+ — iO reside in the fourth 
quadrant while A(e) has a branch cut on [—1; 1]. To pre- 
serve casuality it has to perform an infinitesimal shift up 
where e < and down for e > 0. Such a disposition of 
the poles permits us to choose an integration contour in 
the upper half-plane (the integrand vanishes at infinity as 
e'~ 3 ) and then deform it to a smaller contour C as shown 
in Figure 3. 

Futher, the integrand of Y.fj ial can be split into two parts 
as follows 



u 



de' 



U/2 + iO - A(e') 



U/2 + iO 



One can observe that the second term of the integrand 
does not contribute to the result, since it contains no 




U/2 e + 



Fig. 3: A contour of integration for the constant part of 

singuliarities in the contour. To proceed with the first 
term, we substitute the explicit expressions for ReA(e') 
and ImA(e'). We take into account that only ImA(e') 
changes its sign under a transition from one side of the 
branch cut to another. After routine transformations we 
obtain 



^dual 



'-de' 



U 4t 2 



2 7 rl-4i 2 7_ 1 (e'-e+)(e'-e_) 



(16) 



The integral can be simply done by a trigonometric sub- 
stitution e 1 = sin 6. So we obtain the final answer: 



sr^dual 



U At 2 
2~7r 1 — 4i 2 



7T 

2 



L(e + )-L(e_) 



L(x) = \/\ — x 2 log 



Vl - x 2 + x - 1 
v 7 ! — x 2 — x + 1 



(17) 



(18) 



There is no need to repeat all the calculations for the 
second part (proportional to U 2 ) of equation (11). An 
additional multiplier (e — e' — iO)^ 1 produces another pole 
below the real axis, so it doesn't affect the integration 
contour in any way. This means that the last integral in 
(16) should be replaced with 



de'VT 



(e-e'-i0)(e'-e+)(e'-e_) 
1 f° dt'^Y~^ 



e-e'-iO 



1 



1 



e'-e+ 



In this way we have reduced the J7 2 -part to a known 

'(e) is 

L(e_)-L(e+) , 



result. A full expression for £|| a/ (e) is 



U 



U U 2 { 7T 

2^ 1 - At 2 V 2 "* 
L(e - iO) - L{eJ) 



L(e - iO) - L(e+) 



(19) 



(An infinitesimal imaginary part in Lie — iO) assists to 
choose the right side of the branch cut. Such an imaginary 
addition is not required for e± , since they do not lie on the 
brunch cut for any positive t and U). 

An evaluation of S|" a '(e) repeats the previous one with 
a number of exceptions: 



p-4 



Analytical approximation for SIAM 



• Poles of the integrand are at — U/2 + iO and — e + + iO. 

• The integration contour goes in the opposite direction 
and surrounds [0; 1] segment. 

• The pole of (e' — e — iO)^ 1 lies above the real axis. 
So we have a very similar answer: 

s u a/ ( £ ) = -£rT'(- £ ) (20) 

By taking a limit in formula (19) we achieve even sim- 
plier expression for the t = 1/2 case: 



^fl (ej|t=i/2 — 
eo-U 1 



-idual 



(-e)lt=i/2 = 

L(eo) - L(e - iO) 



l + i(e )-C/ 



1 



eo= 2(^ + 1/^) (21) 

The result for dual self energy should be substituted 
into formula (9). We have plotted thus obtained DOS for 
several values of U and obtained smooth graphs with a 
clear Kondo-like peak at the Fermi level. So, it is clearly 
demonstrated that the simplest first-order theory repro- 
duces the low-energy physics of the Anderson impurity 
model correctly. However, we found that it suffers serious 
problems at higher frequencies. This issue is discussed in 
the next section. 

Renormalization procedure. The mentioned 
drawbacks of the formalism presented so far are related 
to the poles of gw(e). First of all, as it follows from Eq. 
(9), the Green's function G aa is pinned at its mean-field 
value — 1/A at these points. The only case when this pin- 
ning is absent is the vanishing of £^ a ' a * the point of 
the pole, but it seems that no finite-order approximation 
for gw(e) provides such a condition. The pinning is rather 
unphysical, because poles of the atomic problem are by no 
means special points for the entire action (1), and there is 
no any sum rule about them. Further analysis shows that 
the theory fails also near the poles of the atomic prob- 
lem: it replaces the mean- field pole e = U/2 — A(e) with 
the two poles shifted from the real axis in different direc- 
tions. Consequently, the theory is not conservative and 
not causal. 

The problems about the poles of g<jo{e) are probably 
related to the pole structure of the dual Green's func- 
tion. As it follows from the second line of (5), it has two 
close poles, one placed exactly at U/2 and another slightly 
shifted from this point. The residues of these poles have 
opposite signs. We suppose that the description using self- 
energy is not valid in this situation. Anyhow, the above 
paragraph contains a constructive idea on how to improve 
the theory. Since the pinning is absent only if dual self- 
energy vanishes at the pole, 



its reasonable to require the fulfillment of this condi- 
tion. An additional condition means that the theory must 
have an additional ajustable parameter. To introduce it, 
we modify the splitting of the action into Gaussian and 
Hamiltonian parts. We rewrite (1) as follows: 

S = S' at - + fj dtdt'c at (A(t - t') + i\6'(t - t')) c at , 



s 'at = I dt (*(! - A )c ( rt^c CTt - Un n hn + fih at ) 

— OO 

(23) 

The parameter A is to be adjusted to maintain the condi- 
tion (22). 

The calculations are very similar to the above consid- 
ered case A = 0. The atomic problem is the same, up to 
scaling transformations. We obtain 

g aa (e) = ((1 - A)e ± U'/2 T tOy 1 ; U' = (1 - A)" 1 !/ 
j\f n (ei,e 2 ;e 3 ,e 4 ) = -U - C3 X- M 

(24) 

The mean-field consideration gives 

g„(e) = (e±U' -A(e)T^O)" 1 , (25) 

so that the renormalization just results in a shift of the 
atomic resonances. Practical calculation shows that (22) 
is fulfilled with a small positive A, and that A grows as U 
decreases. 

A calculation of the first-order correction is also similar. 
The only important difference is that Q^ al {e) obeys high 
frequency asymptotics yz^e -1 for a finite A, and therefore 
the contour integration contains an additional contribu- 
tion from infinity. The final formula (19) becomes: 

S^ a/ (e;A) = -Eft°'(-e;A) = 

U 4t 2 /_ tt L(e'_) L(e' + ) 
2tt 1 - At 2 V 2 e'_ - e' + 

U' 2 At 2 1 / L(e - iO) - L(e'_ 



2ir 1 - 4i 2 e' + - e'_ 



e — e_ 



L(e-zO)-L(e' + ) \ U> 



G aa (epoic) — 0; g a ^(e poic ) — 0, 



(22) 



_ {U'/2){At 2 -2)±U 2 ^At 2 + {U'/2) 2 -l 
e± - 2(4*2 _ 1} ( 26 ) 

We have performed calculations with this dual self- 
energy substituted in Eq. (9). The value of A has been 
ajusted numerically to fulfill the condition (22). The re- 
sulting graphs for the band DOS are shown in Figure 5. 
Besides band DOS, there are also separated poles, cor- 
responding to the atomic resonances. Their positions are 
shown in Figure 4. For comparison, we also plot the mean- 
field values with and without renormalization (poles of the 
expressions (25) and (5), respectively). It can be verified, 
that the renormalized theory is causal (that is, poles are 
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Fig. 4: Position of the atomic resonance of SI AM with t = D/2 
in different approximations. The resonance bare atom (that is, 
tpoie = tt ) is shown for comparison. Note that the renormal- 
ized theory gives a twice-larger shift of the resonance position 
than the mean-field approach does. For the renormalized the- 
ory, E vanishes near the resonance point, therefore taking the 
dual correction into account almost does not affect the reso- 
nance. Inset shows the value of renormalization parameter A. 



infinitesimally shifted from the real axis in proper quad- 
rants) and conservative (that is, Luttinger theorem is ful- 
filled). 

Conclusion. — At the end of the paper we can com- 
pare our scheme to the approach by Logan et al [9], since 
the local moment approach developed there is formally 
similar to the present theory in a number of points. A 
similar form of the bare Green's functions is used to con- 
struct a diagramatic expansion. But the Green's func- 
tions of the local moment approach are to be determined 
self-consistently, so that it includes poles of a partially 
"dressed" impurity atom. In the same fashion summation 
of ladder diagrams is required to obtain a renormalized 
energy-dependent vertex having a neccesary peculiarity 
near zero frequency. So, the approach by D. Logan re- 
quires an essential numerics. In contrast, we presented 
an almost analytical scheme. Therefore, the result can 
be improved in a regular way by the account of higher- 
order diagrams. Another important peculiarity of our ap- 
proach is that it explicitly describes the renormalization 
of atomic levels due to hybridization with conduction elec- 
trons. A crucial checkpoint for the further development of 
the method is its applicability to multi-orbital systems. 
As far as we know, D. Logan's approach met serious dif- 
ficulties beyond single-band systems, although there is a 
promising result of the variational scheme [14]. For our 
scheme multi-orbital generalization looks straightforward. 

The work was supported by RFFI-DFG grant 08-02- 
91953 — 436-RUS-l 13-938-0. 




Fig. 5: Local density of states p(e) = -SSM ImG(e) of SIAM 
with t — D/2, calculated from the renormalized theory. Atomic 
resonances are shown with vertical lines; height of the line is 
determined by the spectral weight of the resonance. Arrows 
indicate shift from the resonance of the bare atom. 
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